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We have studied the effect of a single nonmagnetic impurity in the recently discovered (K,Tl)Fea;Se2 su- 
perconductors, within both a toy two-band model and a more realistic five-band model. We have found that, 
out of five types of pairing symmetry under consideration, only the dj.2_y2-wave pairing gives rise to impurity 
resonance states. The intra-gap states have energies far away from the Fermi energy. The existence of these 
intra-gap states is robust against the presence or absence of inter-band scattering. However, the inter-band scat- 
tering does tune the relative distribution of local density of states at the resonance states. All these features can 
readily be accessed by STM experiments, and are proposed as a means to test pairing symmetry of the new 
superconductors. 



Introduction. The very recent discovery of high-Tj, (above 
30 K) superconductivity in AFe^Scs (A= K, Tl, Cs) IJll has 
generated a new wave of excitement in the field of iron-based 
superconductors. The new family of compounds is unique in 
several regards, (i) the superconducting transition temperature 
is about four times that of FeSe under ambient pressure (which 
is roughly 8 K) ||4ll; (ii) in contrast to the bad metal behavior 
in other iron-based parent cornpounds, the Fe-deficient com- 
pound (x < 1.6) is insulating yllsn, raising the interest in the 
possibility of Mott insulating state ||6|-(9[| induced by patterned 
Fe- vacancies 1 13, HI; (iii) The end members, TlFe2Se2 and 
KFe2Se2, are heavily electron doped (0.5 electron/Fe) rela- 
tive to other iron-based superconductors (such as LaOFeAs, 
BaFe2Se2, FeSe etc.). Band structure calculations lll2l - ll6ll on 
these end compounds show only electron pockets, primarily 
located around the M point of the Brillouin zone (BZ) defined 
for a simple tetragonal structure. Angle-resolved photoe- 
mission spectroscopy (ARPES) measurements observed these 
electron-like pockets around the M points, and showed no 
hole-like pockets llTl llSll but very weak electron-hke pock- 
ets I 21 ] near the zone center F. 



The superconducting pairing symmetry of iron-based su- 
perconductors has been extensively discussed from both 
weak-coupling and strong-coupling approaches. Within 
weak-coupling approaches, the usual argument for the pop- 
ular s± pairing symmetry relies crucially on the existence 
of F-centered hole pockets and M -centered electron pockets, 
which was the case in previously studied iron-based super- 
conductors. By contrast, strong-coupling approaches invoke 
pairing driven by J1-J2 exchange interactions. The absence 
of F-centered hole pockets, therefore, provides an opportu- 
nity to elucidate the mechanism for iron-based superconduc- 
tivity. Recent calculations have predicted that the supercon- 
ducting state could have dj;2_y2-vt 



s-wave sym- 
metry llgtl, or a mixture of 33.2^,2- and d^j£_^j,2-wave pairing 



states in the magnetic frustration regime Il25ll . All these sce- 
narios lead to nodeless superconducting gap structure, which 
is in agreement with the ARPES observations and other ex- 
periments 271. Because of the particular Fermi surfaces. 



conventional phase-sensitive measurements cannot be readily 
applied to differentiate the different pairing states. 

In this Letter, we propose to use local electronic structure 
around a single nonmagnetic impurity to probe the pairing 
symmetry in (K,Tl)Fei.Se2 superconductors. Such properties 
have proved to be fruitful in identifying the unconventional 
pairing states of different classes of superconductors [12811 . The 
effort has been complicated [29, 2Qi in the previously studied 
Fe-based superconductors due to the existence of both hole- 
like and electron-like Fermi surface pockets. Because the 
Fermi surface of the new (K,Tl)Fej;Se2 compounds comprise 
small pockets of only one type of carriers, this kind of study 
might be very promising. We study the problem within both a 
two-band toy model and an effective model including detailed 
band-structure from LDA calculations. Our results within 
a T-matrix approximation show that the intragap impurity- 
induced bound state only exists when the pairing symmetry 



is of dr 



We also find that the bound-state peak in the 



local density of states (LDOS) occurs at a non-zero energy 
even in the unitary limit, in contrast to the situation in high-Tc 
cuprates. Our prediction can be directly tested by the scanning 
tunneling microscopy experiments on these new compounds. 

Local impurity effects in a toy two-band model. We start 
with a two-dimensional (2D) lattice formed by Fe atoms, and 
consider a two-band toy model with the energy dispersion: 



Em.k = —4:D„i COS /c^i. COS ky + £„ 



(1) 



for each band ?T! (= 1, 2) with the band width 8 1 Dm | and band 
center-of-gravity £,„. In the superconducting state, the bare 
Green's function within the Nambu space is given by 
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where w„ = (2?! + \)'kT (n integer) is the Matsubara fre- 
quency for fermions. The superconducting gap function is 
described by A„i_k- Hereafter, we consider only spin singlet 
pairing. When an alien atom is substituted for Fe, it plays the 




FIG. 1 . (Color online) Energy dispersion along the cut at fcy = (a) 
and the Fermi surface topology (b) in the large one Fe per cell BZ 
based on a toy two-band model with model parameters given in the 
text. 



role of a short-ranged potential scatter In the single-site im- 
purity approximation, the Green's function is given with the 
T-matrix approximation by: 

g{i,j;iujn) ^ go{i,j;iuJn) + Qo{hO;i^^n)TgoiO,j;iuJn) ■ 

Here the T-matrix can be expressed as T{iujn) ~ U[l — 
GaiO, 0; iuJn)U]^^. The potential scattering matrix takes the 
following structure 



U 



(u 





V 







—u 





— V 


V 





u 







— V 





-uj 



(4) 



where u and v are the strength of the intra- and inter-band 
scattering potential. The bare real-space Green's function 
can be obtained from the Fourier transform: QQ{i,j;iujn) ~ 
^o(k; iw„)e*'^'''"'^'"^\ where Nl is the number of 
lattice sites. The LDOS can be evaluated as pi{oj) = 
— ^lm{Ti[Q{i, i;u! + ir])]} with 77 being the intrinsic lifetime 
broadening parameter This quantity is proportional to the lo- 
cal differential tunneling conductance as measured in STM 
experiments. The above scheme is sufficiently general to cap- 
ture essential properties of the single impurity scattering in a 
two-band superconductor 

To mimic the electron pockets of (K,Tl)Fe3;Se2 supercon- 
ductors within the present toy model, we take Di = —1, 
D2 = —0.5, ei = 3, and £2 = —2.5 in the two-band en- 
ergy dispersion as given by Eq. which leads to an electron 
occupation of 2.17 with the zero value of chemical potential. 
Unless specified otherwise, the energy is measured in units of 
\Di\. Figure[T]shows the energy dispersion along the cut with 



and the Fermi surface in the one Fe per unit cell BZ. 



We now turn to the response of the local electronic struc- 
ture to the single impurity scattering in the superconducting 
state with various pairing symmetry. The considered pairing 
symmetry includes conventional momentum-independent s- 



wave channel with A 
channel with Am,k = 
as A„_k = 
nel with A 



- a(") 



Am (sO), the extended s-wave 



cos/s„) (si) as well 



.(SO) ' 

1 


1 






1 

-(s2) 


1 


-rf^i .1 






.(d2) 1 


1- — 1 \ 


1 , \ 


/ , 1 




FIG. 2. (Color online) LDOS measured at the first nearest neighbor- 
ing sites to the single impurity site with only intra-band scattering 
(left column) and with both intra- and inter-band scattering (right 
column). Each column is divided into panels for different types of 
pairing symmetry: (sO), (si), (s2), (dl), and (d2). Except for the 
panel for (dl) in the left column, where u = 1, 5, 100 is taken, 
the scattering strength is taken as u — 100 and v — for the left 
column or u = « = 100 for the right column. The superconducting 



pairing potential is taken as A^'"' - 
broadening parameter 77 = 5 x 10" 



0.08. The intrinsic 



Am-* COS kr,. COS fc„ 



A„j k = Am' sin kx sin ky (d2). For the superconducting gap 
(sO), there is no sign change within each electron pocket or 
across the pockets along the zone boundary. When the radius 
of the electron pockets smaller than tt/2, this is also the case 
for the gap structure (s2). However, the gap structure (si) 
and (d2) have the sign change within each electron pocket. 
The gap structure (dl) has no sign change within each elec- 
tron pocket but changes sign across two neighboring electron 
pockets. The various types of pairing symmetry have been 
considered intensively in the studies of early Fe-arsenic based 
superconductors. 

In our numerical simulation, we have taken the gap func- 
tions in both bands to have same magnitude but out of phase 
with each other We notice that since there are only electron 
pockets from band 1 while the band 2 is below the Fermi en- 
ergy, it is not important whether there is a gap on the second 
band as well as its relative phase with that in the first band. 
Figure |2] shows the LDOS of quasiparticles at a site nearest 
neighboring to the impurity site, which we assume to be lo- 
cated at the center of the square lattice. The left column of 
the figure is for the case without inter-band impurity scatter- 
ing (i.e., V = Q) while the right column is for the case with 
inter-band impurity scattering v = u. The coherent peak 
structure related to the system without impurity substitution 
nicely reflects the momentum dependence of the gap structure 
as just analyzed above. For comparison, we have taken the 



(0) 



(s2), and the d-wave chan- ^^^^ magnitude of the pairing potential aI? for all types of 
(cos/cj: — cosfcj,) (dl) as well as pairing symmetry. Wide quasiparticle gap opens for the (sO), 



(s2) and (dl) pairing states, with their well-like shape being 
similar to what happens in conventional s-wave superconduc- 
tors and compatible to the nodeless gap structure as revealed 
in the ARPES experiments. For the (si) and (d2) pairing sym- 
metry, a y-shape DOS is exhibited reflecting the existence of 
nodal structure in these two types of pairing symmetry. No- 
ticeably, the quasiparticle gap size as measured by the distance 
between two coherent peaks in the DOS is much smaller than 
that for (sO), (s2), and (dl) pairing state. Given the same mag- 
nitude of the pair potential, it has already suggested that the 
superconducting states with (si) and (d2) pairing symmetries 
do not have enough condensation energy; this is consistent 
with the aforementioned microscopic calculations. 

We now turn to the the impurity induced resonance peaks 
when the impurity scattering is in the unitary limit. We find 
that such peaks arise only when the superconducting state has 



the (dl) (i.e., 



-wave) pairing symmetry. To make sure 



that the intra-gap resonance peak originates from the impurity 
scattering, we have also calculated the LDOS by varying the 
impurity scattering potential strength. As shown in the panel 
(dl) of the left column, the resonance peak moves toward the 
coherent gap edge with decreasing potential strength. It is in- 
teresting to notice that even in the unitary limit, the intra-gap 
resonance peak is located far away from the Fermi energy (at 
uj = 0), in contrast to the case in high-Tc cuprates, where 
the imp urity resonance peak is located very close to the Fermi 
energy 112 811 . Two remarks are in order Firstly, the electron 
pockets here are small in size, implying a strong violation of 
particle-hole symmetry when the Fermi energy is located near 
the band bottom. This aspect has been tested to be true for the 
d-wave superconductor as relevant to high-Tc cuprates 128 1. 
Secondly, although there is a sign change of the gap struc- 
ture across two neighboring electron pockets along the zone 
boundary, no nodal quasiparticles (with zero energy) are al- 
lowed in the (dl) pairing state. Inclusion of the interband im- 
purity scattering (see the right column of the figure) does not 
change the above conclusion except that the intensity peak is 
tuned to the opposite side with respect to the Fermi energy for 
a given measure point of the LDOS. 

Local impurity effects in an effective low-energy model with 
realistic band structure. To include the realistic band struc- 
ture at low energies into the study, we have performed band 
structure calculations for KFe2Se2 based on the local density 
approximation. The full-potential linearized augmented plane 
wave (FP-LAPW) method as implemented in the WIEN2K 
code 13111 is used. We then follow the procedure suggest ed by 
Graser et al. 13211 to fit the Wannierized bands 11331 13411 with 
a five-orbital tight-banding Hamiltonian 1351], unfolding the 
small two Fe per unit cell BZ to a large one Fe per unit cell 
BZ. In this procedure, an interface 13611 between the WIEN2k 
code and the wannier90 code 13711 is also employed. As shown 
in Fig. [3] the tight-binding band structure reproduces the LDA 
band structure pretty accurately, and agrees reasonably well 
with those results reported earlier ll22| . Considering the com- 
pound Ko.sFei 7Se2 on which the recent ARPES experiments 
were performed iITtIi . the electron doping is reduced from 0.5 




FIG. 3. (Color online) The five-orbital tight-binding fit (thin solid 
blue lines) of the ten-orbital Wannier fit (thin dashed green lines) 
to the paramagnetic DFT band structure (thick solid red lines). The 
splitting of bands near the Fermi energy at F point has been enlarged 
to guarantee the removal of hole pockets at lower electron doping. 



electron per Fe to 0.1 electron per Fe. Correspondingly, as in 
Ref. |22|, we have artificially adjusted a few tight-binding fit- 
ting parameters to ensure only electron pockets existing in the 
BZ as the chemical potential /i = —0.28 eV is introduced to 
produce a correct filling factor 

The tight-binding Hamiltonian as obtained with the above 
procedure makes it theoretically possible to study the low en- 
ergy quasiparticle properties, where the experimentally ex- 
tracted superconducting gap is at the order of 10 meV. Nu- 
merically, we diagonalize the five-orbital tight-binding model 
with an extremely large number of k points in the large one 
Fe per cell BZ. Out of the 5 bands from the diagonalization, 
the two bands with the one cutting the Fermi energy and the 
other located immediately below the Fermi energy are picked 
out for the calculations of local electronic structure around 
a single nonmagnetic impurity in the superconducting state. 
Guided by the findings from the toy model studied above, we 
consider here only the superconducting state with (dl), that is. 



(i3;2_,y2-wave, pairing symmetry. Figure HI a) shows the den- 
sity of states in the system without the impurity and the LDOS 
at the site nearest neighboring to the impurity site. As shown 
in the toy model, the bare density of states exhibits a well-like 
gap feature around the Fermi energy, indicating a nodeless 
pairing state. Also in the presence of a singe impurity scatter- 
ing in the unitary limit, intragap resonance peak appears in the 
LDOS. The existence of impurity resonance state is also ro- 
bust against the inter-band impurity scattering. In comparison 
to the case of only intra-band impurity potential scattering, the 
location of the resonance energy and its relative intensity on a 
given measure site might be adjusted by the relative strength 
of the inter-band scattering. For a given configuration of im- 
purity scattering, the relative intensity of the resonance state 
depends on the location of the measure site. To better un- 
derstand the profile of the impurity-induced resonance states, 
we have also calculated the spatial dependence of the LDOS 
at the resonance energy, corresponding to the peak locations 
in Fig. lUa). The results are shown in Fig. |5b-c). For the 
case without inter-band scattering, the LDOS intensity is van- 
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FIG. 4. (Color online) LDOS at the 1st nearest neighboring site to the 
single impurity site (a) and the LDOS imaging (b-c) in the supercon- 
ductor with (dl) pairing symmetry. In panel (a), the bare DOS (black 
line) and LDOS without (u — 100 eV and i) = 0) (red line) and with 
(u — V = 100 eV) inter-band scattering are shown. The resonance 
peak Er is at 7.3 meV for the former while 6 meV for the latter when 
the pair potential A^' = 10 meV is chosen. The LDOS imaging is 
taken at oj — Er without (panel (b)) and with (panel (c)) inter-band 
scattering with an intrinsic broadening parameter = 0.5 meV. 



ishingly small at the impurity site regardless of whether the 
selected energy corresponding to the peak at the side above 
or below the Fermi energy (cf. Fig. Hb)). For uj = E^., the 
strongest peak intensity is located at the 1st nearest neighbor- 
ing sites while for uj = —Er (not shown), the strongest peak 
intensity is obtained at the 3rd n.n. sites. In addition, the 
squarish ripples of LDOS distribution is visible for uj = ±Er- 
For the case with the inter-band scattering strength the same 
as the intra-band one, the profile of the LDOS distribution for 
Lo = —Er is very similar to (cf. Fig. Hb)). However, when 
UJ = Er the strongest peak intensity is located at the impu- 
rity site itself (cf. Fig. Stc)) and the LDOS distribution is 
much more localized than the case with the inter-band scat- 
tering. This unique property, associated with the inter-band 
scattering, might enable the usage of the STM technique to 
characterize the relative strength of the inter-band scattering. 

Conclusion. We have studied the effect of a single nonmag- 
netic impurity in the more recently discovered (Tl,K)Fe2Se2 
superconductors within both a toy two-band model and a more 
realistic effective low energy model including band structure 
specifics. We have found that, out of five types of pair- 
ing states under consideration, only the (i^2_j^2-wave pairing 
gives rise to impurity resonance states. The intra-gap states 
have energies far away from the Fermi energy, and they ex- 
ist regardless whether the inter-band scattering is present or 
not. However, the inter-band scattering does tune the relative 
LDOS distribution of the resonance states. Based on these 



results, we propose to use STM measurements of local elec- 
tronic structure near nonmagnetic impurities to probe the p air- 
ing symmetry of the new superconductors. 
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